Quantum Zeno Effect in Radical-Ion-Pair Recombination Reactions 
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Radical-ion pairs are ubiquitous in a wide range of biochemical reactions, ranging from photosyn- 
thesis to magnetic sensitive chemical reactions underlying avian magnetic navigation. We here show 
that the charge recombination of a radical-ion-pair is a continuous quantum measurement process 
that interrogates the spin state of the pair. This naturally leads to the appearance of the quantum 
Zeno effect, explaining a large amount of data on unusually long-lived radical-ion-pairs. 

PACS numbers: 82.30.Cf, 03.65.Yz, 82.20.Xr 



Radical-ion pairs are playing a fundamental role in a 
series of biologically relevant chemical reactions, rang- 
ing from charge transfer initiated reactions in photosyn- 
thetic reaction centers Q to magnetic sensitive reactions 
abounding in the field of spin-chemistry , and in partic- 
ular in the biochemical processes understood to underlie 
the biological magnetic compass of several species having 
the ability to navigate in earth's magnetic field [1,I3|- 

In Fig. [T] we depict a generic model for the radical- 
ion-pair (RIP) creation and recombination dynamics. A 
donor-acceptor molecule DA is photo-excited (D*A) and 
a subsequent charge-transfer creates the RIP (D+A~). 
The singlet and triplet states of the RIP (^D+A", 
•^D+A") are split by internal magnetic interactions of 
the RIP's two unpaired electrons with external mag- 
netic fields and internal hyperfine couplings. The RIP 
is initially created in the singlet state, which is not an 
eigenstate of the magnetic Hamiltonian, and therefore a 
singlet-triplet (S-T) coherent mixing commences. The 
RIP eventually tunnels into an excited state of the neu- 
tral recombined molecule DA, which quickly decays into 
the ground state. As is well known electron trans- 

fer in RIP recombination reactions is fundamentally a 
quantum-mechanical tunneling process. In this Letter 
we will show that this process constitutes a continuous 
quantum measurement of the RIP's spin state. Like ev- 
ery quantum measurement, this one is no exception to 
the rule that measurements performed on a quantum sys- 
tem lead to decoherence However, under appropri- 
ate conditions involving the measurement rate and the 
intrinsic frequency scale of the RIP, the quantum Zeno 
effect Q appears and leads to two physically significant 
consequences: (i) dephasing rates of the RIP S-T coher- 
ent mixing are suppressed, and (ii) the RIP's spin state 
is delocalized, i.e. there is a high probability of triplet 
state occupation even if the singlet and triplet states are 
non-degenerate . 

Quantum Zeno effects appear in several physical sys- 
tems, some of which are very similar to radical-ion-pairs, 
like the ortho-para conversion in molecular spin isomers 
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ultra-cold atom tunneling through optical potentials 
or the suppression of transverse spin-relaxation due 



FIG. 1: Generic level structure and recombination dynamics 
in a radical-ion-pair, taking place in four steps: 1, photoexci- 
tation, 2, RIP creation, 3, tunneling-induced quantum mea- 
surement of RIP's spin state and 4, final decay to the ground 
state. 



to spin-exchange collisions in dense alkali-metal vapors 
[l^ llSj. In the latter case, atomic spin-exchange colli- 
sions, of the form Si • S2, where Si and S2 are the electron 
spins of the two colliding atoms, probe the atomic spin 
state. When the collision rate (measurement rate) ex- 
ceeds the intrinsic frequency scale of the system, which 
is the Larmor frequency of spin precession in the ap- 
plied magnetic field, the effective decay rate of the spin 
coherence is suppressed, a phenomenon that has led to 
the development of new ultra-sensitive atomic magne- 
tometers [IJ]. The RIP tunneling into the neutral state 
is essentially a scattering process [l5|, not unlike atomic 
collisions, that performs a measurement of the RIP's spin 
state, since tunneling can only proceed if the RIP is in the 
singlet spin state. Quantum Zeno effects have been ex- 
tensively analyzed in the literature [11113, H IH , both 
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with respect to pertaining physical systems, as well the 
general conditions leading to the quantum Zeno effect or 
its inverse, the anti-Zeno effect 

In the following, we are going to capitalize on the re- 
markably strong analogy between radical-ion-pairs and 
yet another physical swtem, namely two coupled quan- 
tum dots 23, 24, 25, 2^ 27 1 being continuously interro- 
gated by a point contact. We are going to identify the 
analogous physical observables of the two systems and 
then setup the corresponding evolution equation describ- 
ing the RIP state, in order to directly arrive at the basic 
physical results. The electron hopping between the two 
dots is the analog of the S-T coherent mixing taking place 
in the RIP, whereas the measurement performed by the 
point contact corresponds to the spin-state-dependent 
RIP tunneling into an adjacent excited state of the re- 
combined DA molecule. 

We will consider the simplest possible RIP model, in 
which the triplet-state manifold is degenerate (and de- 
fines the zero energy), and the singlet state has energy 



(1) 



where Q, is the S-T mixing frequency. The tunneling to a 
nearby excited DA state \a) with energy can only oc- 
cur if the RIP is in the singlet state, hence the tunneling 
Hamiltonian is 



'Ht = Tsa cla + Tj 



(2) 



where Tga is the tunneling amplitude. This Hamiltonian 
embodies angular momentum conservation in the tunnel- 
ing process, i.e. the tunneling amplitude for the triplet 
state is zero. In reality, the tunneling Hamiltonian is 
more complicated, since there are several resonant vibra- 
tional states I a), and the recombination rate is given by 
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k = {2nih)\v\^Y.\^sA''K^-^a) 



(3) 



where V is the electronic matrix element and /s,a the 
vibrational overlap between the nuclear wavefunctions of 
a) and the singlet state of the RIP. In this realistic case, 
the tunneling Hamiltonian becomes 



= XI ^Sa 4" + ^Sa ^5°'^ 



(4) 



where now Tsa = l^/s,a- 

Finally, the Hamiltonian of the DA excited state will 
be Tia — Waa^a. The operators cs (cg), ct (c^) and 
a (at) are electron annihilation (creation) operators for 
the single-electron states 15*), \T) and |a), respectively. 
The rate constant k is termed the recombination rate, 
and will be later identified with observable rate con- 
stants. The complete interaction Hamiltonian govern- 
ing the time evolution of the combined system is then 



7i = Ha. + Wrip -I- Ht. The same set of Hamiltonians 
has already been treated at 24|. In similar fashion, it is 



readily shown that by tracing out the \a) degrees of free- 
dom, we arrive at the dissipative evolution of the RIP 
density matrix p: 



dp 
'dt 



-«[Wrip,p] - k'D[clcs]p 



(5) 



where the super-operator acts on the density matrix 
p according to 



V[B]p = B^Bp + pB^B - 2SpSt 



(6) 



The occupation number c^cs can also be written as Qs = 
1/4 — Si -82, which is the singlet-state projection operator 
in a RIP with the unpaired electron spins being Si 2. In 
other words, the eigenvalues of Qs are 1 (RIP in the 
singlet state) and (RIP in the triplet state). Since Qs = 
Qs and — Qs, we arrive at the evolution equation 



dp 
'dt 



-i[HRiP, p] - k[Qs, [Qs, p\] 



(7) 



This is exactly the evolution equation that follows from 
standard quantum measurement theory 0, [1^ , when the 
measured observable is Qs and the measurement rate is 
k. We can generalize this by opening a triplet recombi- 
nation channel, with a recombination rate k^ (the singlet 
recombination rate, so far denoted by k is renamed ks)- 
It is readily shown (since Qs + Qt = 1) that the evolution 
equation is again given by (O, with k = ks + kr- 

In the simple two-dimensional RIP model we are con- 
sidering, the density matrix, the Hamiltonian and the 
singlet-state projection operator are 2x2 matrices: 



\Pts Ptt 



n 



;s = 



1 




(8) 



The evolution equation ([7]) can compactly be written as 
dp/dt — C{p), where £ is a super-operator. The com- 
plex eigenvalues of the matrix A resulting from the re- 
quirement C{p) — Ap (where p is now a column vector 



{pss PsT Pts Ptt) ) ^'^^ of the form —A 



where 



A > is the decay rate and uje is the effective mix- 
ing frequency. In the simple system we are considering, 
there are three independent density matrix elements, and 
hence three eigenvalues. These are shown in Figl2] as a 
function of the recombination rate k (normalized to the 
unperturbed mixing frequency f2). Out of the three de- 
cay rates, two (degenerate ones) increase with k and are 
termed "normal"; the quantum Zeno effect is manifested 
in two ways: (i) the third decay rate, Xqz, is suppressed 
with increasing measurement rate fc, and (ii) the effective 
mixing frequency drops in the limit of strong measure- 
ment (large k). The scaling of Aqz with the measurement 
rate k is Aqz ~ il'^ /k. This scaling is typical of quantum 
Zeno effects , since the Zeno time tz is in this 
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FIG. 2: Decay rates and eigenfrequencies of the evolution 
equation, (a) The "normal" decay rates are the ones that 
increase with the measurement rate k. Quantum Zeno effect 
appears in the decay rate that is suppressed with k, denoted 
by Xqz, which is shown in the inset to be \qz ~ AO? /k at high 
A:.(b) Corresponding eigenfrequencies. The eigenfrequencies at 
high measurement rates are smaller than at low rates, an effect 
well-known from spin-exchange relaxation collisions in dense 
alkali-metal vapors [13] ■ In both (a) and (b) lo/Q, = 1. (c) 
The expectation value (Qs) for p{t = 0) = Qs, lo — O, = 1 
and fc = 0.1 (solid line), k = 1 (dashed line) and = 20 
(dotted line). 



case given by r^^ = V(^I^ripI^) " (^|Hrip|5)2 = f7 
i.e. Xaz = t/t^, where r = 1/k. The long-Uved eigen- 



qZ = t/tI, ' 

value Aqz will determine the long-time behavior of the 
density matrix elements. A typical value of the recom- 
bination rate is fc « 1 ns^^. For an unperturbed mixing 
frequency 11 w 30 ns^^ (this corresponds to a hyperfine 
couphng of 10 G), we find that Aqz ~ 1 /is^^. The hfe- 
time of the S-T mixing process can thus be prolonged by 
several orders of magnitude. 

This naturally explains several experimental observa- 
tions regarding long-lived radical-ion pairs [13, 31 j 
It should be noted that by "recombination" rate we re- 



fer to the fast-decaying eigenmodes of the density ma- 
trix, i.e. those for which the decay rate scales propor- 
tionally with fc. Time scales on the order of 10 ps-1 ns, 
which govern the creation of the RIP through photoexci- 
tation and decay of the excited DA-state could so far not 
be fathomed with long-lived radical-ion-pairs, i.e. since 
the D*A-RIP and RIP-DA energy differences are com- 
parable, why doesn't the RIP disappear at sub-ns time 
scales? This cannot be explained by the presence of the 
"metastable" triplet state, since the mixing rates with the 
singlet state are typically in the 10 ns timescale. There- 
fore, even if the RIP is created in the triplet state, it 
should disappear fast through T-S mixing and singlet- 
channel charge recombination. The quantum Zeno effect 
naturally leads to RIP lifetimes that extend even to the 
/Lts timescale. The theoretical models that were used until 

now p, s, m. Hi] to describe RIP recombination dy- 
namics masked the presence of the quantum Zeno effect, 
since they treated the tunneling process of the RIP and 
the subsequent decay to the DA ground state (steps 3 and 
4 m Fig. [J) with a single, phenomenological density ma- 
trix equation dp/dt = — ^[TYrip, p] — k{pQs + Qsp)- This 
equation accounts for the depopulation of the singlet RIP 
state at a rate 2k, i.e. the probability S — Tr{pQs} to 
find the RIP in the singlet state decays exponentially at 
a rate 2fc, unavoidably inducing a simultaneous dephas- 
ing of the S-T coherent mixing at the the high rate fc. It 
is noted for completeness that the aforementioned semi- 
classical density matrix equation leads to similar results 
just in this particular two-dimensional toy-model of the 
RIP. As soon as we move to a realistic description of 
the RIP, which involves at least an 8-dimensional den- 
sity matrix (4 is the two-electron spin multiplicity and 2 
the spin multiplicity of one nucleus with spin- 1/2, which 
is the bare minimum needed to form a RIP supporting 
singlet-triplet mixing), all decay rates scale proportion- 
ally to fc and there is no manifestation of the Zeno effect. 
This is not the case with the density matrix equation ([7]) . 

The second physical consequence of the quantum Zeno 
effect in the RIP recombination merits some discussion. 
The probability for significant triplet state population is 
high, on the order of unity, even for S-T energy differ- 
ences Lo > n. This effect has been discussed in [23|, and 
in the simple RIP model we are considering, it is seen 
in the fact that pss tends to 1/2 in the long-time limit, 
irrespective of the problem's frequencies. This has tan- 
gible consequences, since in realistic systems, the triplet 
state RIP can recombine to other chemical products, or 
in cases of RIPs in solution, the D and A molecules will 
eventually diffuse away. In both large probabil- 

ity of populating the RIP triplet state will be evident in 
the reduced yield of recombined DA molecules. 

We will finally elaborate on step 4 of FigdJ namely 
the decay to the ground state of the neutral DA molecule. 
The continuous quantum measurement performed by the 
tunneling process into the excited DA state will at times 
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be interrupted when there is a definite measurement out- 
come. This is described by the quantum-jump approach 



24 l37l | of the quantum-trajectories description of dissi- 
pative quantum systems. When the outcome of measur- 
ing Qs is 1, that is, the RIP is in the singlet state for 
sure, tunnehng and decay to the DA ground state can 
proceed. This is formaUy described by the conditional 
evolution of the RIP's quantum state j^"), 

\^{t + dt)) = l'J'(t)) - idt'HRIp\^{t)) 

-kdt{Qs-{Qsmit)) 

+ (9) 

The stochastic point process that takes the values and 
1, i.e. (diV)2 = dN, and M[dN] = 2k{Qs)dt, where M[.] 
represents the mean over all possible realizations of the 
process. The first three term in represent the no-jump 
deterministic evolution of the RIP state, while the last 
term describes the quantum jump that eventually occurs 
opening the possibility for charge recombination, after 
which the DA molecule can relax to its ground state. 
This way we have closed the excitation-recombination 
cycle pictured in Fig. [TJ Analogously to the tunneling 



current [2J, [26| we can define a "charge-recombination 
current" Rc = dN/dt = k{Qs)- The recombination 
reaction rates can then be defined as the characteris- 
tic rates appearing in the two-time correlation function 
G(t) = E[Rc{t)Rc{t + r)], which is similarly to 0, 
shown to be given by G{t) = fc^(Tr{e^'^ps,oo} — 1/4), 
where ps. oc = QsPocQs is the singlet projection of 
the steady state density matrix (which in this case is 
Poo — 1/2). Thus the reaction rates are the eigenvalues 
of C that we have calculated and plotted in Fig. 1. 

In summary, starting from first principles, we have here 
demonstrated the fundamentally quantum-mechanical 
nature of the radical-ion-pair recombination process. 
This provides a natural explanation of several unusual ex- 
perimental findings in RIP chemical reactions. More im- 
portant, it is not inconceivable that quantum Zeno effects 
in radical-ion-pairs could be found to be intimately in- 
volved with the quantum-mechanical foundations of pho- 
tosynthetic reactions. 

I acknowledge helpful discussions with Dr. D. Anglos 
as well as helpful comments by the anonymous Referees. 
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